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When an electric field is applied to an electrolyte-saturated polymer gel embedded with 
charged colloidal particles, the force that must be exerted by the hydrogel on each par- 
ticle reflects a delicate balance of electrical, hydrodynamic and elastic stresses. This 
paper examines the displacement of a single charged spherical inclusion embedded in 
an uncharged hydrogel. We present numerically exact solutions of coupled electrokinetic 
transport and elastic-deformation equations, where the gel is treated as an incompress- 
ible, elastic Brinkman medium. This model problem demonstrates how the displacement 
depends on the particle size and charge, the electrolyte ionic strength, and Young's 
modulus of the polymer skeleton. The numerics are verified, in part, with an analytical 
(boundary-layer) theory valid when the Dcbye length is much smaller than the particle 
radius. Further, we identify a close connection between the displacement when a colloid is 
immobilized in a gel and its velocity when dispersed in a Newtonian electrolyte. Finally, 
we describe an experiment where nanometer-scale displacements might be accurately 
measured using back-focal-plane interferomctry. The purpose of such an experiment is 
to probe physicochemical and rheological characteristics of hydrogel composites, possibly 
during gelation. 



1. Introduction 

Hydrogels are soft, water-saturated networks of polymer with molecular-scale porosity. 
They find widespread use in molecular-separation technologies (gel electro phoresis), drug 



deliv e ry, scaffolds for tissue engineering, cell culture, wound care (e.g.. lLowman et al., 
2004 lYui et"aZll2004T l. and microfluidic pumping and control (e.g.. IBassetti et aZ.ll2005h T 



This paper concerns a novel class of hydrogel composite s where charged col loidal in- 
clusions are immobilized in an uncharged hydrogel matrix. [Matos et al. ( 20061 ) recently 



demonstrated that polyacrylamide hydrogels doped with silica nanoparticles significantly 
enhance the electric-ficld-induced transport of both charged and uncharged molecules 
through the composite. The underlying electroosmotic flow mec hanism is co nsistent with 
theoretical expectations for weak perturbations to equilibrium ( Hilll 2006c ). 



Each charged inclusion experiences electrical, hydrodynamic and mechanical contact 
forces, while the mobile co unter charge produces an electro-osmotic flow that permeates 



the surrounding polymer. IHilll (|2006q ) recently developed an electrokinetic transport 
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model that quantifies how the charge and size of the inclusions, the ionic strength of the 
electrolyte, and gel permeability influence the electro-osmotic pumping capacity of an 
uncharged polymer network with randomly dispersed, impenetrable inclusions. To a first 
approximation, the elastic distortion of the network does not influence the flow or ion 
fluxes, so little is known of the particle displacements and flow-induced distortion of the 
polymer. 

The particle displacement could be used as a diagnostic to probe the physicochemi- 
cal characteristics of the particle-polymer interface, much like the electric-ficld-induced 
particle velocity (electrophoretic mobility) is used to infer the so-called (^-potential of col- 
loidal particles dispersed in Newtonian electrolytes. Also, knowledge of the micro-scale 
strain field in the hydrogel is essential for establishing the electric field strength required 
to initiate micro-scale fracture. Finally, the relationship between particle displacement 
and the elastic and viscous characteristics of the surrounding matrix is central to the 
rapidly advancing field of micro-rheology, which seeks to quantify the dynamics and 
structure of complex fluids ( Solomon fc Lull20"o~Ll Willerbacher fc Oelschlaeger 2007 ). In 
this work, our principal objective is to quantify how the particle, electrolyte and hydrogel 
characteristics influence the electric-field-induced particle displacement. 

As a first step, we solve a model problem where classical electrokinetic transport pro- 
cesses are coupled to the deformation of an incompressible isotropic, homogeneous porous 
medium. The analysis is restricted to situations where the applied electric field is uniform 
and weak, so the particle displacement is small and perturbations to equilibrium may be 
linearized accordingly. In this manner, our approximatio ns are similar to th ose widely 
adopted in the c lassical theories of micro-electrophoresis ( Russel et aZ.lll989h and other 
phoretic motion ( Anderson 19891 ). Nevertheless, our calculations are not restricted by the 
magnitude of the particle ^-potential or the thickness of the equilibrium diffuse layer of 
counterions (Debyc length). 

One promising approach to measure small particle displacements induced by an elec- 
tric field involves collecting light scattered from a micr on-sized inclusion at the focus of 
an optical trap. Back-focal -plane interferometry (e.g., iGittes et aZ.I Il997t ISchnurr et i 



_op 

1997t lAllersma et al\ Il998l ) provides nanometer-resolution position detection with 1 



10 4 Hz temporal resolution. The technique has been successfully applied to measure 
0.1-100 nm displacements of colloidal inclusions in complex fluids (e.g., polymer solu- 
tions and colloidal dispersions) and hydrogels, as well as the electric-field-induced velocity 
(electrophoretic mob ility) of colloi dal spheres disper sed i n Newtonian electrolytes (e.g., 
Galneder et al\ l200l[ ) . Reviews by iMacIntoshl (|l999l ) and iFurstl (|2005l ) present authori- 
tative perspectives of the motivations, development and achievements of micro-rheology, 
with emphasis on single-particle dynamics and the use of optical traps. 

The elastic restoring force of stiff gels significantly attenuates Brownian fluctuations 
in the position of colloidal particulates trapped in the hydrogel skeleton. Accordingly, 
passive micro-rheolo gical techniques are well suited, and indeed limited, to weak gels. 



Schnurr et all (|1997I ) showed that the mean-squared displacement of a particle with radius 



a trapped in an incompressible elastic matrix with Young's modulus £ is kT / (2na£) (kT 
is the thermal energy). With a = 500 nm and £ = 10 kPa, for example, the mean-squared 
displacement is ~ 0.4 nm. Therefore, to register larger displacements in stiffer matrices, 
active micro-rheology typical ly adopts magnetic or o ptical forces to promote oscillatory 
particle motion. For example. I Yamaguchi et al. ( 2005 ) recently demonstrated the use of a 
lock-in amplifier to correlate the particle response in an harmonic optical trap to resolve 
extraordinarily small displacements over a wide frequency range. This proved successful 
for following the transformation (gelation) of viscous polymer solutions to moderately 
stiff visco-elastic hydrogels. 
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Galneder et al. ( 200lh demonstrated the use of electrical forces to facilitate single- 
particle micro-electrophoresis studies. Their experiments produced ~ 1 /im s" 1 particle 
velocities in Newtonian electrolytes at ~ 10-10 2 Hz. Evidently, ~ 10-100 nm particle dis- 
placements were sufficient to measure the electrophoretic mobility. In principle, however, 
smaller electric-field-induccd displacements could be resolved by correlating the response 
(particle displacement) with the applied electric field. 

Accordingly, we complement our principal theoretical results with a description of a 
novel experiment that correlates the forcing and response time series when an electric field 
is used to impart quasi-steady oscillatory motion of a charged sphere embedded in a hy- 
drogel. Because the electric-ficld-induccd displacement reflects the size and charge of the 
inclusions, and the ionic strength of the electrolyte, this experiment could probe physic- 
ochemical changes at the inclusion-hydrogel in terface. These are no t registered using 
conve ntional laser-tweezer micro-rheology {e.g., IValentine et al. 19961 Yamaguchi et al 



20051) because the optical force of the trap, hydrodynamic drag of the electrolyte, and 
the elastic restoring force of an (uncharged) gel are (to a reasonable first approxima- 



tion) i ndependent of the inclusion surface charge (jGittes et 
2001 fef ). 
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The theoretical interpretation of such experiments may be complicated by several 
practical considerations and a variety of interesting physical phenomena. For exam- 
ple, if large electric fields are required to produce measurable particle displacements, 
then non-linear electro-kinet ic influences, such as electroosmosis of the second kind (e.g., 
Mishchuk fc Takhistovll'l995l ) and induced-c harge electroosmosis (e.g.. iBazant fc Squired 
20041 ) come into plav. Aided by experiments. iMishchuk fc Takhistov ( 19951 ) demonstrated 
that electroosmosis of the second kind is significant when the applied (uniform) electric 
field E > 25kT/(ae) (e is the fundamental charge, so kT/e s=s 25 mV at room temper- 
ature). It follows that a linearized electrokinetic model, such as the one pursued in this 
paper, should be reasonable for particles with radius a = 500 nm when E < 13 kV cm" 1 . 
Depending on the electrolyte conductivity, Joule heating may become problematic, so 
lower field strengths are obviously desirable. 

Experiments may also be complicated by electrochemical reactions at the electrodes. 
These influences can be minimized by a dopting time-d ependent (oscillatory) electric 
fields, as exploited in el ectroacoustic (e.g., lO'Brien 1988 ) and dielectric relaxation (e.g., 
DeLacev fc Whitdll98ll ) diagnostic exper iments for colloidal disp ersions, and electrically 
guided colloidal assembly processes (e.g.. iRistenpart et al. 2007 ). Extending the present 
quasi-steady theory to handle oscillatory electric fields at frequencies where dynamics are 
important is a significant challenge. For example, several disparate length scales present 
an extremely 'stiff' computational problem, and dissipativc wave-like dynamics are ex- 
pected from the coupling of viscous (electrolyte) and elastic (polymer) phases. 

Other interesting complications, which are not unique to the propose d electrically 
forced micro-rheolo gy experiments, arise from polymer inhomogeneity (e.g.. lMaggdll998l : 
Huh fc FurstlbOOlU ) and the degree of sticking and slipping at the particle- hydrogel in- 
terface. In principle, these influences can be integrated into the present computational 
methodology with appropriate modifications of the boundary conditions at the particle- 
polymer interface, and by allowing the elastic modulus and hydrodynamic permeability 



of the polymer s keleton to vary with radial position (e.g., iLevine fc Lubenskv! 120016 : 
IChen et al. 2003 ). At present, however, these complications are of secondary importance 
to the task of understanding electrokinetic influences with homogeneous polymer. 

Finally, our analysis concerns uncharged polymer networks. While many hydrogel 
skeletons, especially those of biological origin, bear fixed charge, poly(acrylamide) and 
poly(ethylene oxide) hydrogels, for example, as well as physical gels from block copoly- 
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mers, are, at least in principle, uncharged (jLowman et "aD l2004tlYui et al.ll2004j ). Note 



that electric-field-induced actuation of hydrogels is limited exclusively to poly-electrolyte 
(intrinsically charg ed) hydroge ls, which often respond to relatively weak electric fields 



{E <~ 10 V cm -1 ) ( Shigalll997l ). Although our primary concern is with micro-rheology, 
our theory may contribute to understanding how micro-scale transport processes influ- 
ence macro-scale deflection (electrical actuation) of uncharged hydrogel skeletons that 
are endowed with charge from particulate inclusions. 

Before presenting the full model and examining the results, it is instructive to first 
consider the expected displacement Z of an inclusion if the bare electrical force <jAna 2 E 
is balanced by an elastic restoring force 2na£Z. Here, a is the surface charge density, 
Aira 2 is the surface area, E is the electric-field strength, and £ is Young's modulus of the 
gel. Accordingly, the particle displacement it0 

Z = 2aa£- 1 E. (1.1) 

Equation (|1.1[) over-estimates the displacement by a factor of (2/3)«a when na — > oo 
and |CI < kT/e. Note that na is the particle radius a scaled with the Debye length k^ 1 , 
£ is the (^-potential, kT is the thermal energy, and e is the fundamental charge. When 
na — > and |£| < kT/e, however, Eqn. (jTTTJ) becomes 

Z = 2(e s e £~ 1 E (m -> 0, |C| < kT/e). (1.2) 

Equation (|1.2|) and the correct form of Eqn. (jl.lj) for na — > oo (see section [5]) arc, 
respectively, reminisc ent of the well-know n Huckcl and Smoluchowski limits for the elec- 
trophoretic mobility ( Russel et al. 19891) . Indeed, despite obvious differences, the inclu- 



sion displacement and electrophoretic velocity have a remarkably similar and, in general, 
complicated dependence on the scaled particle size na and scaled ^-potential £e/(fcT). 



2. Models for electrokinetic transport and elastic deformation 

As depicted in figure [TJ we consider a single spherical colloid with radius a and surface 
charge density a embedded in an unbounded, electrically neutral hydrogel. The gel is 
modeled as a homogeneous Brinkman medium that is saturated with an aqueous elec- 
trolyte {e.g., NaCl). Together, the electrolyte concentration, surface charge density and 
particle radius manifest as an electrostatic potential £ at the colloid surface (r = a). Note 
that the counter-charge is concentrated in a diffuse layer with thickness (Debye length) 



2.1. Electrokinetic transport 

The electrokinetic model used in this work to calculate the equilibrium and perturbed 
electrolyte ion concentrations, pressure, a nd fluid velo city is a straightforward exten- 
sion of the standard electrokinetic model (|Booth|[l950l) widely used to describe micro- 
electrophoresis and other electrokinetic phenomena. Details have been presented clsc- 



f Here, the elastic restoring force is the value when a rigid sphere with radius a is embedded 
in an elastic continuum with Young's modulus £ and Poisson's ratio v — 1/2. 
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Figure 1 . Schematic representation of a (negatively) charged, spherical colloid embedded in an 
unbounded electrolyte-saturated, elastic polymer gel (elastic Brinkman medium). An electric 
field E is applied that drives electro-osmotic flow u, which, in turn, exerts a body force on the 
gel that induces a displacement field v. The net result of the electrical, hydrodynamic and elastic 
stress is a displacement Z of the particle (to the left). The dashed lines denote the equilibrium 
(right) and polarized (left) diffuse double layers. In the far-field (r — > oo), the velocity and 
displacement fields decay as r~ 3 . In this work, the displacement field is treated as the sum of 
two fields that each decay as r ; one is induced by a particle displacement Z / with E = 0, 
and the other is the flow-induced distortion with E ^ and Z — 0. 



where ( Hilll 20060 16). The full set of (steady) transport equations is 



e e s V 4> 
D 



N 

£ 



UjZjC 



jU-zjerij — Vip-DjVi 



N 



(2.1) 
(2.2) 
(2.3) 



i]V 2 u - Vp = (i]/£ 2 )u + njZjeVip, 

i=i 

with (steady) ion and electrolyte conservation equations 

V • jj = and V ■ u = 0. 

The electrostatic potential, N ion concentrations and fluxes, electrolyte velocity, and 
pressure are denoted -0, nj and jj, u, and p, respectively. Other variables are the solvent 
dielectric constant e s , pcrmitivity of a vacuum e D , fundamental charge e, ion valances Zj 
and diffusion coefficients £),-, thermal ene rgy fcT, solvent viscosity n, and Darcy perme- 
ability £ 2 (square of the iBrinkmanl (j 1947T ) screening length) . Note that the double-layer 
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thickness (Debye length) is 



K' 1 = ^/kTe s e /(2Ie 2 ) 



where / = (1/2) Ylj=i z j n T * s i° n i c strength, with n°° the bulk ion concen- 

trations. Since this work deals exclusively with steady (or quasi-steady) flows, the fluid 
velocity in Eqn. (|2.3p is relative to a stationary polymer skeleton. 

The equations are solved by perturbing (to linear order) from an equilibrium base state 
where u = and the equilibrium electros tatic potential, ion concentrations and p ressure 
arc denoted ip°, n% and p°, respectively (jO'Brien fc Whitdll978l ; iHill et alj|2003 ). With 



the application of a uniform electric field E, 

ip = ip° + ip', Uj — n° + n'j, u = u' , 

where 

ip' = -E ■ r + $(r)E ■ e r , n'j = hj(r)E ■ e r , 

and 

u = V x V x f{r)E 

= -2] r r-\E ■ e r )e r - (f rr + f r r^ 1 )(E ■ e e )e e . (2.4) 

Here, f r = df/dr, f rr = d 2 f/dr 2 , r is position in a spherical polar coordinate system 
(r, 9, 4>) with unit basis vectors (e r , eg, e^), and e z is the orientation of the polar axis (so 
e r ■ e z = cos 8). 

At the particle-electrolyte interface (r = a), 

u = (no slip), j ■ • e r = (no flux), 

and 

e s e V>^ • e r — e p e V<'!/' • e r = —a (constant surface charge), 

where e p is the particle dielectric constant, and a is the (constant) surface charge density; 
the subscripts attached to the gradient operators distinguish the particle (<) and solvent 
(>) sides of the interface. 
In the far-field (r — > oo), 

ip — ► -E -r + (E ■ e r )D E r- 2 
rij -> nf + (E ■ e r )Cfr- 2 

u -2C E r- 3 (E ■ e r )e r - C E r - 3 (E ■ e 9 )e 9 , (2.5) 

where the scalar coefficients D E , C E and C E n/£ 2 are, respectively, the dipole strengths 
of the electrostatic potential, ion concentrations, and pressure perturbations, induced by 
the electric field. Note that the far-field flow is a Darcy flow, u = — (I 2 /n)Wp' : decaying 
as r -3 as r — > oo. 

The full equations and boundary conditions above are the basis of an efficient numer- 
ical solution that yields the so-called asymptotic coefficients D E , C E and C E , which 
have been used to calcu late the bu lk electrical conductivity and pore mobility of dilute 
polymer-gel composites ( Hill 2006 d fa). This work draws upon C E , and introduces a new 



asymptotic coefficient Z to characterize the far-field decay of the elcctric-field-induccd 
clastic displacement field. 
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2.2. Elastic deformation 

Elastic deformation of the polymer gel is calculated by modeling the polymer skeleton 
as an elastic Brinkman medium, whose equation of static equilibrium is 

V • <x + (i]/e 2 )u = 0. (2.6) 

Here, the elastic stress tensor is (jLandau fc Lifshitd[l986h 

-=(l^[ e+ (l^)( V '^' W 

where v is the (small-amplitude) displacement, £ is Young's modulus, v is Poisson's 
ratio, e = (l/2)[Vt) + (Ve) ], and 8 is the identity tensor. Substituting Eqn. (|2.7[) into 
Eqn. (|2.6|) gives 



^7T^V 2 t;+ f V(V-v) = -(n/f)u. (2.8) 

2(1 + ^) 2(1 + v)(l — 2v) 

Again, the fluid velocity is relative to a stationary polymer skeleton (dv/dt = 0). 

When the (leading order) displacement is divcrgcncc-frcc, which is the situation ad- 
dressed throughout this paper, deformation can only influence the (isotropic) permeabil- 
ity tensor (j]/£ 2 )S by inducing anisotropy. Note that any deformation-induced change in 
permeability yields a non-zero product of the permeability perturbation and the fluid 
velocity, the later of which is itself a perturbation. Accordingly, these second-order terms 
are neglected in the present (linearized) theory. 

There is also a possibility of anisotropy in permeability due to the underlying ran- 
dom microstructure of the medium (hydrogel matrix), which itself is not isotropic at 
very small length scales. However, we assume this medium to be statistically isotropic 
and homogeneous, so that, at the level of a Representative Volume Element (RVE) of 
the deterministic continuum, the an isotropy vanishes just as th e fluct uations in con- 
stitutive response tend to zero; see lOstoia-Starzewski fc Wand (|l989h for a random 



elastic model. Such a scale-dependent homogenization (i.e., a passage from a random 
microst ructure to the RVE) was recen tly studied in the context of Stokesian perme- 
ability ( Du fc Ostoia-Starzewskill2006l). albeit t he de parture from anisotropy was not 
addressed explicitly; see also lOstoia-Starzewskil (|2007j) for related studies in many other 
material problems. 

The Poisson ratios of several widely used, highly swollen, transparent hydrogels are 
reported greater than 0.45. In particular, poly(vinyl alcohol) hyd rogels prepared with a 



mixed solvent of dimethyl sulfoxide and wate r have v ~ 0.472 (|Uravama et al\ 119931 ) 



and polyacrylamide hydrogels have v ks 0.457 ( Takigawa et aLlll996 ). However, it is im- 
portant to note that these measurements arc ascertained from macroscale experiments 
where the characteristic length and time scales cannot probe the equilibrium (long-time) 
state of strain. For example, the fractional change in volume after relaxing to equilib- 
rium is 8V/L 3 r-j (1 — 2v)l/L, where the strain l/L 1 is the ratio of the imposed 
displacement I to the specimen size L. The flux of solvent flowing through the specimen 
during the so-called draining time r is u c ~ 5V/(L 2 t) ~ (1 — 2v)1/t. Further, the flux 
is driven by a pressure gradient p c /L ~ (n/£ 2 )u c , where, to balance the elastic stresses, 
the characteristic pressure is p c ~ £(l/L). Together, the foregoing yield 

t~(1-2v)(i 1 /£)(L/1) 2 , (2.9) 

so with f] ~ 10~ 3 Pas, £ ~ 10 5 Pa, L ~ 10~ 2 m (macro-scale experiment) and I ~ 10 -9 m, 
the relaxation time is r ~ (1 — 2i/)10 6 s. Clearly, this is extraordinarily long if v is not 
sufficiently close to 0.5. Accordingly, when the experimental time scale (e.g., reciprocal 
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frequency) is shorter than the draining time, the change in volume will be smaller than 
at equilibrium, and the apparent Poisson ratio will be greater than the drained value. 

In contrast, the draining time associated with the displacement of a microsphere 
(L = a ~ 10~ 6 m) embedded in a hydrogel with 8 ~ 10 5 Pa and £ ~ 10~ 9 m is 
only r ~ (1 — 2v)10~ 2 s. Clearly, such an experiment is much better suited to prob- 
ing the compressibility of the polymer skeleton. However, solving the problem with a 
compressible polymer matrix demands a distinctly different computational methodology 
than the one adopted in this paper for incompressible gels. Moreover, compressibility 
is anticipated to yield quantitative — not qualitative — changes in the calculated particle 
displacement. This expectation is supported, in part, by the fact that, in the absence 
of an electric field, the quasi-steady particle displacement varies by at m ost 25 percent 
from the incompressible limit as Poisson's ratio spans the range 0-0. (jSchnurr et < 



19971) . Clearly, a definitive answer requires a solution of the electrokinetic equations with 



arbitrary Poisson's ratio: an ambitious goal that is beyond the relatively modest scope 
of the present work. 

Finally, the density and, hence, rigidity of the inclusions (e.g., polymer latex or silica) 
are typically much greater than those of the gel. Accordingly, the inclusions are treated 
as rigid spheres. Moreover, the displacement field is assumed to be continuous across 
the inclusion-hydroge l interface. Other interfacial conditions include the possibility of 
(tangential) slip (e.g. iMura et al. 1985 ) or an opening crack, for example. The opening 
of a crack significantly complicates the analysis, and a slipping boundary condition is 
difficult to justify here, since it requires a physical mechanism to exert an interfacial 
radial stress while maintaining zero (relative) radial displacement and zero tangential 
stress. Accordingly, neither possibility is pursued here. 



3. Superposition to calculate the particle displacement 

It is convenient to calculate the electric-ficld-induced particle displacement Z by su- 
perposing two linearly independent displacement fields. 

One is the displacement field induced by a small displacement Z of the inclusion in the 
absence of an electric field (E = 0). There is no Darcy drag, and the solution of Eqn. (|2.8|) 
can be calculated analytically (see appendix [5}. The resulting mechanical-contact force 
exerted by the polymer on the particle i^j] 

fm ,z = 2ira£Z(l - v) 

1 (5/6-i/)(1 + i/)" 1 ' ' 

The other arises from the Darcy drag force when an electric field E is applied and 
the inclusion is fixed at the origin [Z = 0). The Darcy drag force in Eqn. (|2.8[) is 
calculated from the electrokinetic transport equations with a rigid polymer gel. Then the 
displacement can be obtained by solving Eqn. (|2.8[) . This is detailed in section [4] where it 
is also shown that the mechanical-contact force exerted by the polymer on the inclusion 
is 

f m ' E = [8/3)ttZ e £E - 4ir(ri/£ 2 )C E E (v = 1/2). (3.2) 
In addition to the net mechanical-contact force 

f m = f m ' z + f m ' E , (3.3) 

f Assuming a constant shear modulus £/[2(l + u)] (|Landau~fc Lifshitz 19861 ). 
J With a slipping boundary condition (i.e., zero radial displacement and zero tangential stress 
at r = a), the force is -12n£aZ(l - v)/[(7 - 8i/)(l + v)]. 
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there are electrical and hydrodynamic (drag) forces acting on the particle, denoted f e ' E 
and f ' E , respectively. These are already known from earlier solut ions of the electroki- 



nctic transport equations with a rigid (unperturbed) polymer gel ( Hill 20066l fc). Their 



sum can be written in terms of the asymptotic coefficient C E that characterizes the 
far-field decay of the electric-field-induced flow: 

f e ' B + f d ' E = ATr(r]/f)C E E. (3.4) 

Finally, static equilibrium of the particle demands 

f e,E + f d,E + f m,Z + f m,E = q (35) 

Therefore, collecting the explicit expressions for the various forces above [Eqns. ()3. 1|) . 
p.2p (with v = f/2) and ()3.4j) ] gives the particle displacement 

Z= {4/3){Z E /a)E (1/= 1/2). (3.6) 

The task of calculating Z E is detailed in the next section. Note that the contributions 
involving C E vanish, indicating that the slowest (r^ 1 ) far- field decay of the displacement 
field vanishes upon superposition. In other words, there is no net force acting on the 
polymer, so the Darcy drag force exerted by the electrolyte on the polymer is counter- 
balanced by the mechanical-contact force exerted by the inclusion on the polymer. In a 
composite with a finite inclusion number density, or finite volume fraction </>, part of the 
net mechanical-contact force acting on the polymer must be provide d by a mec hanical 
support to balance an accompanying 0((f>) average pressure gradient (lHillll2006ch . 



4. The electric-field-induced mechanical- contact force for an 
incompressible, elastic Brinkman medium 

This section addresses the displacement induced by Darcy drag when an electric field 
is applied and the inclusion is fixed at the origin. This problem is adopted to calculate 
the force f m - E appearing in Eqn. (|3.2[) . Note that the numerical solution is limited to 
incompressible (V • v = 0) displacement fields, so v = 1/2. 

In appendix[Bj the displacement field is expanded as a power series in a small parameter 
e=l — 2v, i.e., 

v = v Q + evi + .... (4.1) 

The leading contribution to the displacement vq is divergence-free (V • Vq = 0) and 
satisfies the 0(e) equation of static equilibrium, 

V 2 « + V(V-«i) = -(^ 2 )«| (f = l/2). (4.2) 

Incompressibility (as required by the 0(1) problem) is guaranteed by writing 

v = V x V x g{r)E, (4.3) 

where g(r) is a function of radial position r. It follows that 

t> = -2g r r- l {E ■ e r )e r - (g rr + g r r~ l ){E ■ e e )e 9l (4.4) 

where, for example, g r = dg/dr. 

Because the fluid is also incompressible, its velocity field may be written as 

u = V x V x f(r)E, (4.5) 



where f(r) is available from earlier work e xamining th e influence of an electric field (|Hilll 
and a bulk concentration gradient 1200661 ) with a rigid polymer gel (v = 0). 
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If the deformation is assumed not to affect electrokinetic transport processes, which 
is a reasonable approximation when the displacement is divergence-free (so the polymer 
segment density and, hence, the Darcy permeability are unperturbed), then the earlier 
calculations also provide an exact solution with (weak) elastic deformation. 

Substituting Eqns. (|4.3|) and l|4.5|) into the curl of Eqn. (|4.2p gives 



where 



^ S+ §^f = 0, (4.6) 

and f(r) is known. In this manner, Vq is decoupled from V\. 

The fourth-order ordinary differential equation for g r (r) is solved numerically as two 
second-order differential equation^: 

d 2 h 4dh 4 37, d 2 f 2df r 2 

T2 + ~a 2 h = -TFit-rr + 2^> ( 4 - 8 ) 

dr z r dr r z ti dr r dr r z 

ff-fc (4.9) 

When the displacement is continuous across the inclusion-hydrogcl interface, and van- 
ishes in the far-field, the boundary conditions are 

u = at r = o, (4.10) 

and 

v -> -2Z E r~ 1 (E ■ e r )e r - Z E r- 1 (E ■ e e )e e as r -> oo. (4.11) 

The asymptotic coefficient Z E characterizes the strength of the r _1 decay of v (reflecting 
a net force). It follows that 

g r = g rr = at r = a, (4-12) 

and 

g r — ► Z E and g rr — > as r — > cxj. (4-13) 

An analytical boundary-layer analysis that solves the problem when na 3> 1, f < a 
and Id < kT/e serves to verify the numerical solution and to highlight the parametric 
scaling of Z. The result is presented in section [5] [Eqn. (|5.3|) ]. where we also examine 
numerically exact solutions of the full model. 

Turning to the force, the leading contribution to the isotropic stress requires knowledge 
of the 0(e) displacement field Vi, which is not divergence-free (see appcndix[B]). Clearly, 
the divergence of V\ is necessary to evaluate the leading contribution to the force. Fortu- 
nately, the isotropic stress can be obtained from the solution of the 0(e) problem above. 
This is achieved by integrating Eqn. (|B 4[) once i?o is known. The task is simplified even 
further, because only the far-field decay of the displacement and fluid velocity fields are 
needed to evaluate the force. We have verified our general procedure by applying it to 
two simpler problems: one is the classical problem of Stokes flow past a sphere, and the 
other is the elastic restoring force on a rigid sphere embedded in an incompressible elastic 
continuum. Recall, the exact solution of the latter problem, for any v [see Eqn. (|3.1[) ]. is 
worked out in appendix [XJ 

f A finite-difference scheme based on Hill et al.'s methodology (jHill et a/.ll2003f ) is adopted. 
This features an adaptive, non-uniform grid to handle the disparate length scales. 
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The mechanical-contact force exerted by the polymer on the inclusion is 

f m ' E = I crndA = I cr-hdA+ [ (r)/l 2 )udV. (4.14) 

J r—a J r — >oo J r—a 

Note that (see appendix |B|) 

If 1 

<r = Y le + -(V-v 1 )6], (4.15) 



where e = (1/2) [Vv + (Vt> ) T ], from Eqn. ([42]) 

V • vi = - [ [(S7 2 v Q ) + (3/£ ){il/t 2 )u] ■ e r dr, (4.16) 

J oo 

and 

u = -2f r r- 1 (E ■ e r )e r - (f rr + / r r" 1 )( J E • e e )e e 

-> -2C E r~ 3 (E ■ e r )e r - C E r~ 3 {E ■ e e )e e as r ^ oo. (4.17) 

Recall, C E is the asymptotic constant that represents the dipole strength of the far-field 
pressure field (decaying as r~ 2 ) that drives the far-field Darcy flow (decaying as r -3 ). 

Evaluating the first integral on the right-hand side of Eqn. (14. 14|) over the surface of a 
large concentric sphere gives 

I (T-hdA=(8/3)TrZ E £E-(4/3)Tr(i 1 /l 2 )C E E (is = l/2). (4.18) 

J r — >oo 

The volume integral can be transformed to another integral over the surface of a large 
concentric sphere (V • u = and u{r = a) = 0) giving 

(n/£ 2 )udV = {i]/i 2 ) / (x ■ u)hdA = -{8/3)tt(?i/£ 2 )C e E. (4.19) 

-a J r—>oo 

Finally, adding Eqns. (|4~T8]) and |4~T9)l gives the mechanical-contact force as it appears 
is Eqn. (^2j) . 



5. Results 

When solving the equations numerically, the characteristic scales adopted for length, 
velocity and displacement are 

k~ , u* = e s e (kT/e) 2 /(r]a) and rju*/£ = e s e (kT/e) 2 /(£ a), 

respectively. It is therefore convenient to introduce a dimensionless asymptotic coefficient 
Z E so 

e (kT/e )a , 

£( Ka y 

Accordingly, the particle displacement [Eqn. (|3.6[) ] is 



ZE = c)l 2 



e.«e (kT/e) ~ E 



The independent dimensionless parameters adopted below are na, C,e/(kT) and k£. 
Note that Z E is independent of £ , and, furthermore, we will see that the displacemen t 
is a very weak function of £. Using a boundary-layer approximation for C (|Hillll2006a[ ). 
an analytical solution for Z E , when ra > 1, f <C a and \(\ < kT/e, yields 

Z — > 3(e s e £~ 1 E as na — » oo, (5-3) 



12 R. J. Hill and M. Ostoja-Starzewski 

which is the counterpart to Eqn. (|1.2|) identified in the introduction. 



5.1. Particle displacement with NaCl electrolytes 

To draw a closer connection to experiments, Z/E is plotted in figures [5] and [3J with 
£ = 1 kPa. Therefore, actual displacements Z (nm) can be conveniently obtained by 
multiplying the ordinates by the electric field strength E (V cm -1 ) and dividing by 
Young's modulus £ (kPa). 

Figure shows Z/E as a function of the scaled ^-potential £e/(kT) for a particle 
with radius a = 500 nm and a hydrogel with Young's modulus £ — 1 kPa (v = 1/2) 
and Brinkman screening length £ = 5 nm. The electrolyte is NaCl, with ionic strengths 
corresponding to na — 1-10 3 . As expected, the (negative) particle displacement is in 
the direction of the electrical force, 4ira 2 aE, and increases with the magnitude of the 
C-potential (C < 0). 

At low ^-potentials, the displacement is clearly proportional to \(\ and, hence, the sur- 
face charge density a = e s e nQ (when \Q\ < kT/e). As expected, the numerical solutions 
approach the boundary-layer theory [Eqn. (|5.3|) ] as na — > oo. However, the numerical 
results reveal distinct maximums at moderate and large values of |£|; these are due 
to polarization of the d iffuse double layer. As is well known from electrophoresis {e.g., 
O'Brien fc White] 1978 ). polarization diminishes the local electric field, thereby atten- 



uating the electrical force. Because polarization by electro-migration and relaxatio n by 



molecu lar diffusion are practically independent of the polymer gel in this model (jHill 
2006 d ). they are as significant here as they are in electrophoresis. 



Figure [3J shows Z/ E under the same conditions as in figure [H but now as a function 
of the scaled reciprocal double-layer thickness kcl, with each curve corresponding to a 
constant ^-potential. The displacement is clearly a weak function of na, particularly 
when | C| is small, and, as expected from figure [2 a much stronger function of £. Clearly, 
this way of presenting the results emphasizes the large values of na required for the 
boundary-layer theory to be accurate. 

Since figures [5] and [3] are presented with a fixed value of i/a = 0.01, it remains to 
establish the influence of the Darcy permeability (or Brinkman screening length). Recall, 
the boundary-layer theory [Eqn. (|5.3|) ] indicates that the displacement is independent of 
I. More generally, however, the particle displacement reflects the far-field electric-ficld- 
induccd distortion of the polymer skeleton when the particle is fixed at the origin by an 
external force 

f E = -f dE - f eE - f m,E = -{8/3)itZ e £E. (5.4) 

Figure |4] shows how the particle displacement varies over four decades of the scaled 
Brinkman screening length k£ when £e/(kT) = —6 and na = 10~ 2 -10 3 . While there are 
obvious transitions when k£ ~ 1 (from plateaus where k£ — > and oo) , the displacement 
is remarkably insensitive to I. 

In the next section, we identify a simple transformation that permits Z E and, hence, 
the particle displacement, to be approximated by the well-known electrophoretic mobility, 
which, of course, is independent of I. Accordingly, we write 

Z = 3Ce s e £ _1 /[™, Ce/{kT), k£]E, (5.5) 

where the dimcnsionless function 



(5.6) 
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Figure 2. The ratio of the displacement Z (nm) to the electric field strength E (V cm ') as a 
function of the scaled ^-potential C,e/{kT) for various scaled reciprocal double- layer thicknesses 
ko = 1, 10, 100 and 1000: NaCl at T = 298 K; a = 500 nm; I = 5 nm; £ = 1 kPa; v = 1/2. 
Note that the displacement is inversely proportional to Young's modulus £. The dashed-line is 
the boundary- layer theory [Eqn. (|5.3[) ]. 



and, to a reasonable approximation, 

f[Ka,(e/{kT),Kl] w f[Ka,(e/(kT),oc] 



(5.7) 



It follows that figure [5] or [3] is sufficient to span a significant range of the parameter 
space (strictly for negatively charged inclusions and NaCl electrolyte). Furthermore, from 
Eqns. HH]) and (|S~3]) . it is evident that / -> 2/3 as na -> (with \C\< kT/e) and / ->• 1 
as na — > oo. 



5.2. Connection to the electrophoretic mobility 

The displ acements shown in figure s[2]and[31bear a close resemblance to the electrophoretic 
mobility (jO'Brien fc Whitelll978t) . so we present in figures [5] and [5] a scaled displacement 



(Z/E)£e/{2e s e kT) = (3/2)/[«a, (e/(kT), Ki](e/(kT) 



(5.8) 



for a symmetrical electrolyte (KC1) (solid lines, figure [5]) and a representative asymmetric 
electrolyte (Ba(NOa)2) (figure [S]). This dimensionless quantity (solid lines) has a very 
similar dependence on na and Qe/{kT) as the scaled electrophoretic mobility (dashed 
lines) 

(U/E)3 V e/(2e s e kT) (5.9) 
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Figure 3. The ratio of the displacement Z (nm) to the electric field strength E (V cm 1 ) 
as a function of the scaled reciprocal double-layer thickness na for various scaled (^"-potentials 
-Ce/(fcT) = 1, 2, 4, 6 and 8: NaCl at T = 298 K; a = 500 nm; I = 5 nm; £ = 1 kPa; z/ = 1/2. 
Note that the displacement is inversely proportional to Young's modulus £ . The dashed- lines 
are the boundary- layer theory [Eqn. (|5.3p ]. 
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Figure 4. The ratio of the displacement Z (nm) to the electric field strength E (V cm 1 ) as 
a function of the scaled Brinkman screening length kI with (e/(kT) = —6 and ko, = 0.01, 0.1, 
1,..., 1000: KC1 at T = 298 K; a = 500 nm; 8 = 1 kPa; ^ = 1/2. Note that the displacement is 
inversely proportional to Young's modulus £. 
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Figure 5. The scaled displacement — (Z/E)£e/(2e a e kT) as a function of the scaled ("-potential 
—C,e/{kT) for various scaled reciprocal double-layer thicknesses na = 0.01, 0.1, 0.2, 0.4, 1 and 
2 (solid lines, left panel); and ko = 3, 6, 10, 15, 30, 50, 100 and 1000 (solid lines, right panel): 
KC1 at T — 298 K; a = 500 nm; £ — 5 nm; v = 1/2. Dashed lines are the scaled electrophoretic 
mobility (U / E)3rje/(2e a e a kT) for selected values (labels in parentheses) of na — 0.01, 0.1, 0.2, 
0.4, 1 and 2 (left panel); and 10, 30 and 100 (right panel). 



presented in O'Brien and White's well-known paper (see O'Brien fc White1l978l figures 3, 
4, & 6jJ Here, U is the electrophoretic velocity, i.e., the translational velocity acquired by 
a charged spherical colloid dispersed in a Newtonian electrolyte in response to a (weak) 
electric field E. 

The difference between the scaled mobility and scaled displacement highlighted in 
figure [S] is small when \(\ < kT/e, but increases appreciably when na ~ 1 and |£| > 
3kT/e. Note that the scaled mobility is smaller than the respective scaled displacement, 
because the electrophoresis problem over-estimates the convective contribution to the ion 
fluxes, thereby over-polarizing the diffuse double layer and, therefore, under-estimating 
the electrical force in the particle-displacement problem. 

A simple but approximate relationship between Z/E and U/E can be established by 
eliminating the Darcy drag force from the fluid momentum equation [Eqn. (|2.3p ] and 
the polymer equation of static equilibrium [Eqn. (|2.8[) ]. This produces the momentum 
conservation equation in the electrophoresis problem [i.e., Eqn. (|2.3[) without the Darcy 



f The electrophoretic mobilities reported here were calculated using software (called MPEK, 
ava ilable from the corresponding author) based on the methodology of Hill, Saville and Rus- 
sel (jHill et aLl2003l ) for the electrophoretic mobility and other single-particle properties of spher- 
ical polymer-coated colloids. Here, the influence of a polymer coating is removed by specifying 
an infinite (large) Darcy permeability for the coating. The m obilities in figures [5] an d [I?] are 
in excellent agreement with O'Brien and White's calculations (jO'Brien fc White] 1 19781 ); small 
differences can be attributed to the opposite sign of the (-potential. 
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Figure 6. The scaled displacement —(Z/E)£e/(2e 3 e kT) as a function of the scaled (^-potential 
—C,e/(kT) for scaled reciprocal double-layer thickness kb = 5: KCl and Ba(NOa)2 at T = 298 K; 
a = 500 nm; £ = 5 nm; ^ = 1/2. 



drag term] with a modified 'fluid' velocity 

ii = u + v£/(3rj). 



(5.10) 



Recall, u is the fluid velocity in the polymer gel and v is the displacement of the skeleton. 

By setting u = u' in the ion conservation equation, the solution of the electrophoresis 
problem (involving only it') over-estimates the convective ion fluxes by an 0[(v c /u c )£ / (3r))] 
amount; here, v c and u c are, respectively, characteristic polymer displacement and fluid 
velocity scales. However, since the convective ion fluxes are O(Pcj), where the Peclet 
numbers Pcj = u c k~ 1 jDj 1, the absolute errors are small. 

Note that the polymer displacement reflects a transfer of the electrical body force from 
the fluid to the elastic skeleton. This transfer occurs by direct coupling of the fluid and 
polymer (via the Darcy drag force) and through indirect coupling by the transfer of vis- 
cous stresses from the fluid to the particle, which, in turn, are transferred to the polymer 
via mechanical contact between the particle and polymer. Consequently, the polymer dis- 
tortion must be independent of the permeability in so far as the electrical body force is 
constant. However, the electrically drive n flow incre ases significantly with the permeabil- 
ity [either as I or £ 2 , depending on k~ - 1 (lHillll2006clfah ]. so the 0[(v c /u c )£/(3rf)] errors m 
the ion conservation equations must diminish with increasing permeability. Accordingly, 
the solution of the electrophoresis problem must yield an increasingly accurate solution 
of the particle displacement problem as k£ — > oo. 

More quantitatively, the solution of the (E) electrophoresis problem ( O'Brien fc White] 

19781) yields u' ~ C Er 1 as r — > oo. Therefore, because v ~ Z E Er 1 and u ~ C E Er 3 
as r — > oo, it follows from Eqn. (|5.10p that 



C E 3r]/£ as n£ 



oo 



[y = 1/2). 



(5.11) 
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Furthermore, since the electrophoretic mobility 

M = U/E = C E /C U 7 (5.12) 

— E — U 

where C and C are the asymptotic coefficients associated w ith the far-field decay o f 



the fluid velocity in the (E) and (U) (electrophoresis) problems (jO'Brien &: Whitdll978f ). 
Eqns. (EU), (joTTTj) and (|5\T2l) give 

Z F(3r)/£ )ME as nl ->■ oo (f=l/2). (5.13) 

Note that F ~ (4/3)C U /a = D s /D ~ 1 is the particle drag coefficient, which can be 
conveniently expressed as the ratio of the Stokes- Einstein- Sutherland diffusivitjlj D s = 
kT/(6Trna) to the actual diffusivity D[na,(e/(kT)] > D s . The drag coefficient is shown 
in figure [7] for negatively charged colloidal spheres in a KC1 electrolyte. Since 1 < F < 1.2 
and most often 1 < F < 1.02, the error in the particle displacement that comes from 
setting F = 1 in Eqn. (|5 . 13[) tends to be small compared to the error arising from finite 
k£. 

To demonstrate the correctness of Eqn. (|5.13[) . let us briefly consider a specific example 
where £e/ (kT) = —6 and na = 10. From figurc[4l the scaled displacements are ~ 3.89 and 
« 3.42 when k£ = 0.1 and n£ — * oo, respectively. Furthermore, from figure [5l the scaled 
mobility and drag coefficient are ~ 3.36 and « 1.016, respectively. Therefore, multiplying 
the scaled mobility by the drag coefficient gives w 3.41, which, as expected, compares 
extremely well with the value (« 3.42) obtained directly when k£ 3> 1. 

In summary, we have established that the electrophoretic mobility and drag coefficient 
combine [as indicated in Eqn. (|5.13[) ] to yield the correct limiting value of the scaled 
particle displacement as k£ — > oo. Furthermore, the difference between the displacement 
with impenetrable and infinitely permeable polymer reflects the degree to which poly- 
mer influences convective ion transport. Finally, because polarization by convection is 
weak relative to electromigrative polarization and diffusive relaxation, the influence of 
permeability on the displacement is generally small. It should be noted, however, that 
Young's modulus (or shear modulus) of the skeleton and the permeability both vary 
with the polymer density according to scaling laws t hat have been studied extensively in 
the polymer physics literature (e.o.. lde Gennej 1979h . In practice, therefore, particle dis- 



placements are expected to vary primarily with the modulus, with relatively insignificant 
changes due to the accompanying change in permeability. 



6. Application to micro-rheology 

As demonstrated above, dimensional particle displacements arc < 10e s e o kTE/ (e£) 
in the experimentally accessible range of na and \Q\e/(kT). In polymer gels saturated 
with aqueous electrolytes at room temperature, this implies that Z < 0.l7(E/£) nm. 
For example, with £ = 10 kPa and E = 100 V cm -1 , we find Z < 0.17 nm, which 
poses a significant, but not necessarily insurmountable, experimental challenge. When 
£ = 0.1 kPa, however, Z < 1.7 nm is achieved with E = 10 V cm -1 . 

As highlighted in the introduction, passive micro-rhe ology also tends to be limited to 



comp lex fluids whose storage modulus is < 1 kPa (e.g.. ISchnurr et ali 119971 : IChen et al 



20031 for F-actin solutions and polyacrylamide gels, and A-DNA, respec tively). For hydro 



gels, such conditions are realized when the cross-linking density is low (|de Genneslll979h . 



As our theory demonstrates, high electric field strengths (> 1 kV cm : ) and/or low 

t ISauires fc Bradvl (|2005f l present a compelling case to associate W. Sutherland with this 
famous relationship. 
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Figure 7. The drag coefficient F = (4/3)C /a = D s /D[na, C,e/{kT)] for charged spherical 
colloids as a function of the scaled (^-potential —C,e/{kT) for various scaled reciprocal double-layer 
thicknesses na = 0.01, 0.1, 0.2, 0.4, 1 and 2 (left panel); and na = 3, 6, 10, 15, 30, 50, 100 and 
1000 (right panel): KC1 at T = 298 K. 



elastic moduli (< 1 kPa) are necessary to register nanometer-scale electric-field-induced 
particle displacements. Accordingly, to minimize the electric field strength and, hence, 
avoid non-linear electrokinctic phenomena and Joule heating, we propose an experimental 
methodology that seeks to resolve sub-nanometer particle displacements. 

In contrast to existing laser tweezers micro-rheology, where t wo lasers are used to 
simultaneously perturb and image a micron-sized tracer particle Valentine et fl L » : 



Yamaguchi et a/1l2005l ) , we envision an experiment that necessitates only an imaging laser 



to measure Z/E = (4/3)(Z E /a), which, recall, reflects physicochemical characteristics 
of the particle-gel interface. In practice, an optical trap and piezo-electric stage (with 
feed-back control) may be adopted for positioning and alignment purposes. 

Very small-amplitude particle displacements induced by an oscillatory electric field may 
be registered by correlating the measured response Z' = Z— (Z) = Z"(t) + \Z\ cos (ut + 0) 
with the (harmonic) electric field strength E = \E\ cos (tot). Here, Z"(t) denotes the noise 
(adjusted to have zero mean), \Z\ and \E\ are the amplitudes of the (harmonic) particle 
displacement and electric field, 6 is the phase angle, u> is the angular frequency, and t is 
time. If the noise Z" is uncorrelated with the excitation E, then 

(Z'E)^^ 1 [ Z'(t)E(t)dt -> (l/2)|Z||^|cos6» + 0(T- 1 / 2 ) (6.1) 



if the duration T of the time series is large. 

The quasi-st eady particle equation of motion [Eqn. (13.511 with viscous drag force 
—6irr]*adZ/dt ( Valentine et al. 19961 : Yamaguchi et al. 20051) ] yields a (complex) dis- 
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placement 



with amplitude and phase given by 

\Z\ te/gg^ and (6-3) 

v / i+P?W 

Therefore, if the viscosity and elastic modulus arc known, from a complementary macro- 
or micro-rhcology experiment, then 

(4/3)(Z £ /a) = — f'E) 

(1/2) cos6»|£'| 2 v /l + (tanfl) 2 

Moreover, at sufficiently low frequencies (i.e., cu <C £/(3r)*)), the quasi-steady response 
is simply 

(4/3){Z E /a) M 2(Z'E)/\E\ 2 (v = 1/2). (6.5) 

In this manner, it is not necessary to have explicit knowledge of the visco-elastic charac- 
teristics of the gel to directly probe the physicochcmical characteristics of the particle, 
as characterized by (4/3)(Z E /a). Furthermore, as a compelling alternative to a more di- 
rect measurement of Z/E to determine (A/3)(Z E / a), the correlation (Z'E) may provide 
an effective means of dealing with the small signal-to-noise ratio expected in such an 
experiment. 

Finally, when the frequency is higher than the reciprocal draining time [Eqn. (|2.9p j. a 
non-trivial extension of the quasi-steady theory is necessary to account for a compress- 
ible polymer skeleton and o scillatory relative motion of the fluid and polymer phases. 



Levine k, Lubenskvl (|2001d ) undertook such a study for particles embedded in polymer 
gels without electro kinetic influences. They sign ificantly extended the prevailing quasi- 
steady theory (e.g., Levine fc Lubenskvl 12001 ij ). identifying the frequencies where, in 



particular, compressibility of the polymer skeleton influences the dynamical response. 
Clearly, their work provides a sound basis for advancing the present theory to elucidate 
dynamical electrokinetic effects in micro-rhcology applications. 



7. Summary 

We presented a theoretical model to calculate the electric-field-induced displacement of 
a charged, spherical colloid embedded in an electrolyte-saturated polymer gel. The stan- 
dard electrokinetic model describes electric-ficld-induced transport of ions and electrolyte 
momentum, with a Darcy-drag term that couples the electrolyte mass and momentum 
conservation equations (Brinkman's equations) to a continuum equation of static equi- 
librium for an unbounded, incompressible, linearly clastic polymer skeleton. 

The scaled particle displacement 

(Z/E)£e/(2e s e kT) 

has a similar dependence on fe/(fcT) and na as the well-known scaled electrophoretic 
mobility 

{U/E)3r]e{2e s e kT). 

More precisely, we showed that the product of the scaled electrophoretic mobility and 
particle drag (friction) coefficient yields the exact scaled displacement when k£ — > oo 



20 



R. J. Hill and M. Ostoja-Starzewski 



(i.e., when the polymer skeleton presents zero hydrodynamic resistance to flow). How- 
ever, because the particle displacement decreases only slightly as n£ passes through ~ 1, 
when increasing from to oo, the scaled electrophoretic mobility provides an excel- 
lent approximation of the scaled displacement over a wide range of the experimentally 
accessible parameter space. Therefore, as expected from electrophoresis, the scaled dis- 
placement is linear in the ^-potential when |£| < 2fcT/e, and has distinct maximums 
when |e| ~ 5kT/e. The decrease in displacement with increasing ^-potential is due to po- 
larization and relaxation of the diffuse double layer. Evidently, polarization is dominated 
by elcctromigration, since 1 the influence of convection, which can be completely arrested 
by an hydrodynamically impenetrable polymer skeleton, is extremely weak. 

To help deal with the small signal-to-noise ratio expected in experiments that ne- 
cessitate sub-nanometer particle displacements, we established a connection between 
Z/E = (4/3)(Z E /a), which reflects the physicochemical state of the particle-gel interface, 
and the correlation of the response (Z') and forcing (E) signals in a back-focal-planc in- 
terferometry apparatus. As a particle characterization tool, the novel micro-electrokinetic 
experiment we describe is analogous to classical micro-clcctrophorcsis, with existing laser 
tweezer micro-rheology serving as the (complementary) analogue of dynamic light scat- 
tering. 

We did not examine the stresses in the surrounding polymer. Nevertheless, our cal- 
culations provide an important first step toward future studies aimed at quantifying 
the micro-scale states of stress (and strain) when these soft composite materials are sub- 
jected to electric fields. Our model may be helpful for studying fracture, and it provides a 
means of quantitatively interpreting experiments designed to measure small electric-ficld- 
induced displacements of charged inclusions. In turn, these could be used to probe the 
mechanics of weak (uncharged) polymer gels at length and time scales that are beyond 
the reach of conventional (macroscalc) rhcomctcrs. 

RJH gratefully acknowledges support from the Natural Sciences and Engineering Re- 
search Council of Canada (NSERC) (grant number 204542) and the Canada Research 
Chairs program (Tier II). 
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Appendix A. The force to displace a finite sized sphere embedded in 
a compressible elastic continuum 

This appendix provides an analytical solution of the equation of static equilibrium 
in the absence of body forces (Darcy drag). In turn, the force required to displace (by 
distance Z) a rigid sphere (with radius a) embedded in an unbounded elastic continuum 
(with Young's modulus £ and Poisson's ratio v) is obtained. 

Substituting a solution of the form 

v = v +vi, (Al) 

where 

V 2 v = 0, (A 2) 

into the equation of static equilibrium [Eqn. (|2.8[) ] gives 

VV + _ X V[V ■ (v + m)] = 0. (A3) 

Taking the curl yields 

V 2 (Vxv 1 ) = 0, (A 4) 

which, with the prevailing axial symmetry, provides a scalar equation for (V x Vi) ■ e.$. 
Since V x V\ is an harmonic pseudo vector, symmetry and linearity yield a non-zero 
decaying solution V x v\ ~ Z x Vr _1 . It follows that v\ ~ Zr _1 , which is harmonic 
and, hence, can be attributed to vq. Accordingly, the only non- harmonic contribution to 
Vi is irrotational (V x v\ = 0) and, hence, 

vi=V<j>, (A 5) 

where <fi is a scalar function of position. Substituting this into Eqn. (|A 3j) gives 

V[V • v + 2(1 - j/)V 2 0] = 0, (A 6) 

so 

V ■v a + 2(l-is)V 2 (j) = 0. (A 7) 

Again, symmetry and linearity considerations yield the general decaying solution 

v = c Z + aZr' 1 + c 2 (Z ■ V)V(r- x ). (A 8) 

It follows that 

V • v = Cl Z ■ V(r _1 ) (= -dZ ■ e r r- 2 ) (A 9) 
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and, hence, 

Note that V 2 ip = r 1 has the solution ip = r/2, so writing Eqn. (|A 10p as 

V ^=- 2(1-,) (AU) 



requires 



and, hence, 



ciZ • e r 

'4(1-1/) 



(A 12) 



01 = = - 4 Cl _ [Z - (Z • e r )e r ]. (A 13) 

Finally, the complete displacement field is 

w = cqZ+ ^Z - - Cl [Z - (Z ■ e r )e r ) + %(Z • e r )e r - Z], (A 14) 
r 4r(l — v) r 3 

where the scalar constants Co, Ci and C2 must be chosen to satisfy the boundary condi- 
tions. 

A.l. No-slip boundary condition 
If, for example, v — > — Z as r — > oo and t> = at r = a (fixed), then cq = —1 and 

= Z + °±Z + ^§[3(Z • e r )e r - Z] - - Cl \Z - (Z • e r )e r ], (A 15) 
a a J 4a(l — v) 

which requires 

6(1 — v)a o? , , _ „. 

Cl = lb6iy andC2 = 2(6^- (A16) 

The mechanical-contact force on the inclusion is therefore 

Jr=a (1+1/) (5-6l/)(l+!/) 1 ' 

Note that the displacement field can be rewritten in terms of the force, so the Green's 
function 

G= J^[(3-4,) 5 + e r e r ] (A 18) 

is obtained by changing reference frames (v = Z at r = a, and t; — » as r — > oo) and 
letting a — > 0. 

A. 2. S7ip boundary condition 

Again, if v — > Z as r — > oo, but u ■ e r = at r = a (zero radial displacement) and 
er • e r — (er • e r ) • e r e r = at r = a (zero tangential stress), then cq = —1, and 

6(1 -v)a , (l-2,)a 3 /A 

ci = '— and c 2 = '—. A 19) 

(7-8,) 2(7-8,) v ; 

The mechanical-contact force (on the inclusion) is then 

2^£ciZ _ 12n£aZ(l-u) 
T ll + v) ~ (7-8,)(l + ,)' (A2U) 
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Appendix B. The leading-order isotropic stress for an incompressible 
elastic continuum 

Writing the displacement field as 

v = vo + ev\ + e 2 V2 + (Bl) 
where e = (1 — 2v) <C 1, substituting this into the equation of static equilibrium 

V 2 v + -V(V ■ v) = -(rj/£ 2 )u^^-, (B2) 

and collecting terms of like order in e gives at 0(1): 

V(V-u ) = 0; (B3) 



W 2 v + V(V-vi) = -fa/^u-; (B4) 



VV+V(V ■v 2 ) = {j 1 /t 2 )u\; (B 5) 



at 0(e) 
at 0{e 2 ) 
and at 0(e 3 ) 

V 2 v 2 + V(V • v 3 ) = 0. (B6) 

Note that, if u and Vq are both divergence- free, then V • Vi, V ■ v 2 , etc., all satisfy 
Laplace's equation with general solution, e.g., V • v\ = a\ + b\E ■ Vr _1 . 

Because the leading contribution to the displacement v is divergence-free [Eqn. (jB 3jl ]. 
it can be written 

» = VxVxj(r)f;. (B7) 
Substituting this into the curl of Eqn. jB4j gives 



where 



= < B8 > 

v 2 = ^S» < B9 > 

and /(r) is known. Clearly, the solution is independent of Vi. Note, however, that V\ 
contributes to the leading-order stress, 

If 1 

a = —[e + -(V-v 1 )8}, (BIO) 

where e = (1/2) [V«o + (Vt)o) T ]. It follows that the leading contribution to the integral 
of the surface traction is 

f'o = \ [ [V^o + (Vv f + (V • • e r dA, (B 11) 

where [Eqn. (|B4|) ] 

V • t;i = - / [fo/* 2 )«4 + V 2 -y ] • e r dr'. (B 12) 

Note that only the far-field decays of u and Vo are necessary to evaluate this integral 
when r — > oo. Recall, 

it -> -2C E r~ 3 (E ■ e r )e r - C E r~ 3 (E ■ e e )e e as r ^ oo, (B 13) 
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so 



V • «i -> [2Z E - (r)/l 2 )-C E ]r- 2 (E ■ e r ) as r -> oo 



(B14) 



and, hence, 



/0= o 



3 7, 



e r jdA 



= £(16/9)nZ E E + £(8/9)nZ E E - (4/3)ir(r)/£ 2 )C E E 

= (8/3)irZ E £E - {4/3)w(ri/£ 2 )C E E (B 15) 

Finally, adding the volume integral [-{8/3)ir(r]/£ 2 )C E E] gives the net mechanical-contact 
force acting on the inclusion. 



errt,E 



{8/3)ttZ e £E - 4ir(r]/t)C E E 



(B16) 



